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Abstract 



We report on the bright Lyman-break galaxies (LBGs) selected in a 767 arcmin 2 area of the Subaru 
Deep Field. The selection is made in the i — zr versus zb — zr plane, where zb and zr are new bandpasses 
with a central wavelength of 8842A and 984lA, respectively. This set of bandpasses enables us to separate 
£ — \ well z ~ 6 LBGs from foreground galaxies and Galactic cool stars. We detect 12 LBG candidates down 

CO 1 to Z K = 25.4, and calculate the normalization of the rest-frame far-ultraviolet (FUV: ~ 1400A) luminosity 

function at Aipuv = -21.6 to be 0(-21.6) = (2.6±0.7) x 1CT 5 mag" 1 Mpc~ 3 . This must be the most reliable 
measurement ever obtained of the number density of bright z ~ 6 LBGs, because it is more robust against 
both contamination and cosmic variance than previous values. The FUV luminosity density contributed 
from LBGs brighter than M F uv = -21.3 is (2.8 ± 0.8) x 10 24 ergs s _1 Hz" 1 Mpc" 3 , which is equivalent 
■ to a star formation rate density of (3.5 ± 1.0) x 10~ 4 A/ Q yr _1 Mpc~ 3 . Combining our measurement with 

those at z < 6 in the literature, we find that the FUV luminosity density of bright galaxies increases by 
q ■ an order of magnitude from z ~ 6 to ~ 3 and then drops by 10 3 from z ~ 3 to the present epoch, while 

the evolution of the total luminosity density is much milder. The evolutionary behavior of bright LBGs 
resembles that of luminous dusty star- forming galaxies and bright QSOs. The redshift of z ~ 3 appears to 
be a remarkable era in the cosmic history when massive galaxies were being intensively formed. 

Key words: galaxies: evolution — galaxies: high-rcdshift — galaxies: luminosity function, mass 
function — galaxies: photometry 



1. Introduction Very recently, several groups have made surveys of a 

similar galaxy population at z ~ 6 using ultra-deep imag- 

Over the last decade, selection of star-forming galaxies ing data taken with Hubble Space Telescope (HST) and 

at high redshift using continuum breaks seen in their far- 8m-class telescopes. Most of them search for galaxies red 

ultraviolet (FUV) spectra has been applied successfully in i — z. Such 'i-dropout' galaxies are likely star-forming 

to deep multi-color images taken with large telescopes, galaxies at z^ 6, since the continuum break of z^ 6 galax- 

Thousands of galaxies beyond z ~ 2 have been detected ies falls into the i band. Systematic studies of galaxies at 

by this method in two-color planes; galaxies selected in z ~ 6 and beyond are crucial not only for the undcrstand- 

this way are called Lyman-break galaxies (LBGs). For ing of the formation and early evolution of galaxies during 

instance, LBGs at z ~ 4 can be isolated from foreground and just after the reionization epoch of the universe, but 

objects (lower-redshift galaxies and Galactic stars) in the also for the identification of objects which contributed to 

B — R vs R—i plane, since they have extremely red B — R the reionization. 

colors and relatively blue R — i colors. Various properties We begin with a brief summary of the previous studies 
of LBGs up to z~ 5 have been discussed in detail based on on z ~ 6 galaxies. Fontana et al. (2003) made a multi- 
photometric data and follow-up spectroscopy (e.g., Steidel color Very Large Telescope (VLT) and HST imaging sur- 
et al. 1999; Papovich et al. 2001; Shapley et al. 2001, vey aiming at detecting 4.5 < z < 6 galaxies. In their VLT 
2003; Giavalisco & Dickinson 2001; Porciani & Giavalisco data for an area of about 30 arcmin 2 , 13 candidates with 
2002; Foucaud et al. 2003; Iwata et al. 2003; Ando et al. Zab — 25 mag were found based on photometric redshifts, 
2004; Ouchi et al. 2004a,b). among which four were at 5 < z < 6. Dickinson et al. 
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(2004) detected in early Advanced Camera for Surveys 
(ACS) GOODS data of 316 arcmin 2 five z ~ 6 candidates 
with i— z > 1.3 down to z ~ 26 mag. Stanway et al. (2003, 
2004) detected eight and six i-dropout (i — z> 1.5) galax- 
ies with z < 25.6 in public HST ACS data of the GOODS- 
South (146 arcmin 2 area) and the GOODS-North (200 
arcmin 2 ), respectively. Bouwens et al. (2003) detected 
23 z ~ 6 galaxy candidates in two ACS Guaranteed Time 
Observation fields (46 arcmin 2 in total) down to z ~ 27.3 
mag. A deeper survey was made subsequently by Bouwens 
et al. (2004a) for two Hubble Ultra Deep parallel fields; 
30 i-dropouts (i — z > 1.4) were found over 21 arcmin 
down to z = 28.1. Bouwens et al. (2004b) discussed the 
size evolution of galaxies up to z ~ 6 in the Hubble Ultra 
Deep Field (HUDF) and the Hubble Ultra Deep parallel 
fields. Yan, Windhorst, & Cohen (2003) found 30 z ~ 6 
candidates brighter than z~ 28 mag in an ACS/WFC par- 
allel field (10 arcmin 2 ). Yan & Windhorst (2004) found 
108 candidates with z < 30.0 mag in the HUDF (10.34 
arcmin 2 ). Bunker et al. (2004) also made a search in 
the HUDF, finding 54 candidates at z < 28.5 mag with 
i — z > 1.3. Many of the studies have reported a decrease 
in the FUV luminosity density of bright galaxies by a 
factor of w 3 — 8 from z ~ 3 to z ~ 6 (Dickinson et al. 
2004; Giavalisco et al. 2004; Stanway et al. 2003, 2004; 
Bouwens et al. 2004a; Bunker et al. 2004), although 
Fontana et al. (2003) and Bouwens et al. (2003) obtained 
a roughly constant luminosity density over this redshift 
range. Measuring the number density of bright LBGs ac- 
curately is of great importance in discussing and modeling 
the formation and evolution of massive galaxies. 

There are, however, two main uncertainties involved in 
the previous studies of bright z ~ 6 LBGs. First, selection 
of LBGs by a single i — z color is likely to suffer from a 
significant amount of contamination, especially at bright 
magnitudes, since red galaxies at z ~ 1 — 2 and Galactic 
cool stars can be as red in i — z as z ~ 6 LBGs. Removal of 
unresolved objects in a high-resolution HST image is an 
efficient way to reduce the contamination by cool stars. 
However, high-resolution imaging is not helpful to sepa- 
rate red z ~ 1 - 2 galaxies from LBGs. Near infrared (NIR) 
photometry is useful especially for discriminating LBGs 
from red galaxies. In the previous studies, however, NIR 
photometry was either limited to a very small sky area, 
or not deep enough to securely separate red galaxies from 
LBGs. Second, the areas surveyed in the previous z ~ 6 
LBG papers are a few tens to a few hundreds arcminutes, 
the widest being 346 arcmin 2 . This may not be sufficiently 
large when one recalls an extremely low surface density of 
bright LBGs (an order of ~ 10 -2 arcmin" 2 ) and large cos- 
mic variance in their volume density expected from obser- 
vations of lower-redshift LBGs and biased galaxy forma- 
tion models. Thus, wide-field and deep two-color imaging 
is essential for constructing a large and clean sample of 
z ~ 6 LBGs. 

In this paper, we report on a survey of z ~ 6 LBGs 
which is aimed to overcome the two main uncertainties in 
the previous surveys. We apply for the first time a two- 
color selection to imaging data as wide as 767 arcmin 2 . 



We obtain a reliable value of the number density of bright 
z ~ 6 LBGs. By combining our value with those in the 
literature, we examine the evolution of the number density 
of bright LBGs over 3 < z < 6 and discuss the formation 
of massive galaxies. 

The plan of this paper is as follows. Section 2 explains 
the strategy of our survey. Section 3 describes the photo- 
metric data used in our analysis. Selection of z ~ 6 LBGs 
is made in §4. Section 5 is devoted to results and dis- 
cussion, and a conclusion is given in §6. AB magnitudes 
are used in this paper. We assume a flat universe with 
Q = 0.3, A = 0.7, and H = 70 km s" 1 Mpc" 1 . 

2. Strategy of Our Survey 

Our survey covered an area of 767 arcmin 2 of the Subaru 
Deep Field (SDF: Maihara et al. 2001). This wide survey 
area will considerably smooth out large-scale cosmic vari- 
ance which will be probably at present at z ~ 6. Several 
observations have reported large cosmic variance for z > 5 
galaxies. Bremer et al. (2004) found that the surface 
density of z ~ 5 galaxies varies considerably over areas 
comparable to a single ACS pointing. Shimasaku et al. 
(2003, 2004) reported large inhomogencity of the sky dis- 
tribution of Lyman a emitting galaxies at z ~ 5 on scales 
of 50 Mpc. At z ~ 6, a transverse comoving length of 50 
Mpc corresponds to 20'. Hu et al. (2004) and Ouchi et 
al. (2005) also found a highly inhomogeneous spatial dis- 
tribution of ~ 5.7 Lyman a emitting galaxies on tens of 
Mpc scales. See Somerville et al. (2004) for a theoretical 
prediction for cosmic variance of high-redshift galaxies. 

The bandpasses we use for two-color selection are i, zb, 
and zr. The zb and zr filters, which were designed by us 
for this purpose, divide the SDSS z band into two at 9500 
A. Figure 1 plots the effective response curves of the zb 
and zr bandpasses which include the filter transmission, 
CCD sensitivity, throughput of the optics of the Subaru 
Telescope, and atmospheric extinction at an airmass of 
1.2. The central wavelength and FWHM of the zb and 
z R bandpasses are (8842 A, 689 A) and (9841 A, 537A), re- 
spectively. These two bandpasses enable us to measure the 
slope of continuum emissions redward of Lyman a, which 
is crucial for discriminating z ~ 6 objects from z ~ 1 - 
2 red galaxies and cool stars both of which have similar 
i — zr but redder zb — zr colors (see §4). 

We have another three advantages of using zb and zr. 
First is that we can measure reliably rest-frame UV con- 
tinuum flux from zr photometry The continuum flux of 
i — z selected LBG candidates is usually measured from 
their z-band magnitude. However, a recent study based 
on spectroscopic observation of an z-dropout galaxy at 
z = 6.33 suggests that a large fraction of i-dropout galax- 
ies may have a strong Lya emission (Nagao et al. 2004). 
If this is true, measurements of the continuum flux of 
z > 6 galaxies from z-band photometry would be biased 
high due to contamination from Lya emission. A galaxy's 
Lyman a line does not enter the zr bandpass as long as 
the galaxy is at z < 6.7. 

Second, we can construct an LBG sample of a well- 
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Fig. 1. Response functions of the zb and zr bandpasses 
(solid lines). The two dashed lines indicate the response of 
the i and z bands, respectively, of the Subaru/Suprime-Cam. 

defined and relatively narrow redshift range thanks to 
zb — color. As will be shown in §4. the zb — zr value 
of star-forming galaxies increases rapidly with redshift at 
z > 6. Thus, selection of LBG candidates with zb — zr 
color less than an appropriate value (zb — zr = 0.7 in this 
study) provides an LBG sample of a redshift distribution 
with a sharp upper limit, although the precise correspon- 
dence between zb — zr and redshift varies with the intrin- 
sic spectral shape of LBGs. This is important in calcu- 
lating the comoving volume surveyed by our observation. 
In contrast, the selection function of a single i — z color 
criterion is wide, having a long high-redshift tail. 

Third, our z ~ 6 LBG sample is free from contamination 
by spurious objects near the detection limit of the zr data. 
As will be seen in §4, one of our selection criteria for z ~ 
6 LBGs is zb — zr < 0.7. Since the zb data are deeper 
than the zr data by about 0.7 mag, this criterion ensures 
that objects satisfying this criterion are detected in the zr 
image as well. Fake objects in the zr image which happen 
to exceed the detection flux are basically undetected in the 
Zr image, and thus will be eliminated by this criterion. 
Selection by i — z color alone could be contaminated by 
spurious objects near the detection limit. 



Table 1. Summary of the observational data. 
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777.1* 
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26.2 


2 July 2003 
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25.4 
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public 
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27.7 


public 
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600 


27.8 


public 
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801 


27.4 


public 
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504 


26.6 


public 


NB816 


600 


26.6 


public 


NB921 


899 


26.5 


public 


* 3<7 limiting magnitudes on 


a 2" 


iiamctcr aperture. 



3. Data 

3.1. zb and zb. Imaging 

We took deep zb- and ZR-band imaging data of 
a 30' x 24' area in the Subaru Deep Field (SDF) 
[13 h 24 m 38. s 9, +27°29'26"(J2000)] with the prime focus 
camera (Suprime-Cam; Miyazaki ct al. 2002) mounted on 
the 8.2m Subaru telescope on 2003 July 2 and 2004 March 
19 and 20. The image scale of Suprime-Cam is 0."202 per 
pixel. The exposure time was 85 min for zb and 599 min 
for zr; we set a much longer exposure time for zr, since 
the effective throughput of the zr band is considerably 
lower than that of zb due to the rapidly decreasing quan- 
tum efficiency of the CCDs toward longer wavelengths be- 
yond 9500 A. The individual CCD data were reduced and 
combined by IRAF and a data reduction software pack- 
age for Suprime-Cam developed by Yagi et al. (2002) and 
Ouchi (2003). The combined images for individual bands 
were aligned and smoothed with Gaussian kernels so that 
the PSF sizes of the final images match those of the pub- 
lic images described in the next subsection. The final 
images after removal of low S/N regions near the edges 
of the field of view cover a contiguous 767 arcmin 2 area 
with a PSF FWHM of 0."98. Photometric calibration was 
made using the spectrophotometric standard star GDI 53 
(Bohlin, Colina, & Finley 1995) taken on 2003 July 2 for 
both bandpasses. A summary of the data is given in Table 
1. 

3.2. Public Imaging Data 

The SDF has deep, public Suprime-Cam data of seven 
bandpasses, B, V, R, i, z, NB816, and NB921, obtained 
in the Subaru Deep Field Project (Kashikawa et al. 2004). 
This project is a large program of Subaru Observatory to 
carry out a deep galaxy survey in a blank field in optical 
and NIR wavelengths to study distant galaxies. B, V, 
i?, i, and z are standard Johnson and SDSS broad-bands, 
and NB816 and NB921 are narrow-bands whose central 
wavelength and FWHM are (8150 A, 120 A) and (9196 
A, 132 A), respectively. The exposure times and limiting 
magnitudes of the public images are given in Table 1. We 
combine the i-band data with our zb and zr data to select 
LBGs at z ~ 6 in the i — zr vs zb — zr plane. We also use 
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Fig. 2. Number counts of galaxies in the zjjband. Detection 
incompleteness has not been corrected. The error bars imply 
Poisson errors. 

the B,V,R data to further remove foreground objects from 
the LBG sample. 

3.3. Object Detection and Photometry 

Object detection and photometry were made using 
SExtractor version 2.1.6 (Bertin and Arnouts 1996) on all 
nine images (zb, zr, and seven public images). The zr- 
band image was chosen to detect objects. If more than 
5 pixels whose counts were above the 2 <7 s k y were con- 
nected, they were regarded as an object. In total, 45,405 
objects were detected down to zr = 25.4 mag (3a limiting 
magnitude). Figure 2 plots the number counts of objects 
detected in the zr band, where 1 stellar objects with the 
PSF shape brighter than zr = 23 have been removed. For 
each object detected in the zr image, a magnitude on a 
2"-diameter aperture was measured for each passband to 
derive the colors of the object; we adopted for each object 
the same aperture among all passbands to ensure that the 
colors were measured correctly. We adopted MAG_AUTO 
for the total zr magnitude. 

4. Selection of z ~ 6 Lyman-break Galaxies 

Figure 3 shows the selection of LBGs at z ~ 6 in the 
i — zr vs zb — zr plane. The magenta curves in this figure 
indicate the tracks of young, star-forming galaxies at z > 5 
with three dust extinction values, E(B — V) =0,0.15,0.50; 
the filled squares on the curves correspond to z = 6. The 




Fig. 3. Distribution of objects in the i — 2p_ versus 23 — zr 
plane. The dots indicate objects detected in our data. If the 
i or 2b magnitude of an object is fainter than the 2cr mag- 
nitude of the corresponding bandpass, it is replaced with the 
2<r magnitude. The dots along the line i — 2r = 23 — 2r + 1.25 
at 2b — zr > 0.8 are objects which are fainter than the 2a 
magnitude in both i and 2b- The 12 LBG candidates are 
shown by filled circles (i magnitudes brighter than the 2cr 
mag) and filled triangles (fainter than the 2cr mag). The ma- 
genta dashed lines represent the predicted tracks of young 
star-forming galaxies with three different dust extinction val- 
ues: from lower left to upper right, E(B — V) = 0,0.15,0.5. On 
the tracks, redshifts are marked by small squares with 0.2 in- 
terval; the filled squares correspond to 2 = 6. The green lines 
imply the tracks of elliptical galaxies redshiftcd from 2 = to 3 
without evolution; the spectra are taken from Fukugita et al. 
(1995; solid line) and Coleman et al. (1980; dashed line). The 
cyan crosses represent Galactic stars given in Gunn & Strykcr 
(1983) and Knapp et al. (2004). The adopted boundary for 
LBG selection is shown by the thick solid line. 

value of E(B — V) = 0.15 is the median extinction of z ~ 4 
LBGs detected in the SDF (Ouchi et al. 2004a). For the 
intrinsic spectrum of the young galaxies, we adopt the 
spectral synthesis model of Kodama & Arimoto (1997) 
and set an age of 0.1 Gyr, Salpeter IFM, and a star- 
formation timescale of 5 Gyr. We use Calzetti et al.'s 
(2000) formula for dust extinction and Madau's (1995) 
prescription for the absorption by the IGM. It is found 
that the i — zr colors of the model galaxies increase mono- 
tonically with rcdshift up to z = 6 while their zb — zr colors 
stay nearly constant. At z > 6, the model galaxies become 
rapidly redder in zb — zr. 

The crosses represent 175 Galactic stars of a wide va- 
riety of types given in the spectrophotometric atlas of 
Gunn & Stryker (1983) and 46 Galactic cool stars taken 
from Knapp et al. (2004) for which spectra are available 
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down to 7000 A. The green solid and green dashed curves 
represent the tracks of elliptical galaxies from z = to 
z = 3 (without evolution) given in Fukugita, Shimasaku, 
& Ichikawa (1995) and Coleman, Wu, & Weedman (1980), 
respectively; these two curves are plotted to show the po- 
sition and its typical uncertainty in the i — zr vs zb — z r 
plane of the reddest interlopers. Coleman et al.'s spec- 
trum does not show a UV upturn and thus has redder 
colors than Fukugita et al.'s. 

Figure 3 demonstrates that LBGs at z ~ 5.6 - 6.2 
are separated well from foreground objects (intcrmcdiate- 
rcdshift red galaxies and Galactic stars) in the i — zr vs 
zb — plane. A simple selection of LBGs using a single 
color, e.g., i — zr, may suffer from contamination, since 
red galaxies at z ~ 1 - 2 and cool stars have similar i — zr 
colors to z ~ 6 LBGs; zb — zr, color plays an essential role 
in discriminating such interlopers from LBGs. 

Previous two-color selections of z 5 LBGs have usu- 
ally adopted for the vertical axis of the two-color plane 
(i.e., the axis sensitive to the Lyman break) the bluest 
bandpass minus the middle-wavelength bandpass; for ex- 
ample, for S,i?,i-selected LBGs at z^4, B — R, not B — i, 
is adopted as the vertical axis. In this study, however, 
we have adopted i — zr, not i — zb, for the vertical axis 
from the technical reason that since we are based on a 
ZR-detected catalog, some objects are undetected in zb, 
for which i — zb colors cannot be expressed properly. Of 
course, adopting i — zb instead of i — zr does not change 
our result as long as we use the same selection criteria. 

From Figure 3, we define the selection criteria for LBGs 
at z ~ 6 as: 

i Zr > Zb Zr + 1.3, 

zb - z R < 0.7, 
«-zr>4(z b -z r ) + 0.1, 
B,V,R>2a mag. 

The selection boundary in the i — zr vs zb — zr plane de- 
fined by these criteria is shown by the black solid lines 
in the figure. These criteria select a model galaxy with 
E(B - V) = 0.15 at z = 5.6 - 6.2. Varying E(B-V) over 
and 0.5 does not significantly change the redshift range se- 
lected. To further reduce contamination from foreground 
galaxies, we impose the fourth criterion that an LBG can- 
didate must be un-detected, i.e., less than the 2a limiting 
flux, in any of the B, V, and R images. 

The dots indicate the objects detected in our data. 
The black circles and triangles indicate LBG candidates 
in our data which satisfy all the selection criteria; trian- 
gles correspond to objects not detected in i, i.e., fainter 
than the 2a magnitude. All the observed magnitudes and 
colors have been corrected for Galactic absorption using 
Schlegel, Finkbeiner, & Davis (1998). 

We select 12 LBG candidates in total down to zr = 25.4, 
which are summarized in Table 2. We do not plot the error 
bars in Figure 3 to avoid crowding. The typical la errors 
in zb — zr and i — zr for the 12 candidates are 0.3 mag 
and 0.4 mag, respectively (The i— zr error is estimated for 
objects with i magnitudes brighter than the 2a mag). The 



total z R magnitudes (MAG_AUTO) of all 12 candidates 
are fainter than 24.8. However, all the candidates are 
visible in both zb, z, and NB921 images above the 3a 
level, meaning that they are not spurious objects. 

The FWHMs of the 12 candidates in the Zr image span 
over 0."86 - l."78, with a median of l."2 and an interquar- 
tile range of l."14 - l."53. Similar FWHM values are ob- 
tained for the z image which is considerably deeper than 
zr. For most of them, the deviation from the PSF is not 
significant, when we consider that SExtr actor tends to re- 
turn FWHM values slightly larger than the PSF size for 
very faint point sources. Visual inspection of the 12 can- 
didates in the z image finds that three objects seem to be 
resolved, with a major axis of about 2". The relatively 
large seeing size of our data prevents us from measuring 
the intrinsic sizes of candidates, but some brightest LBGs 
at z ~ 6 could have a large star forming region of several 
kpc. 

Snapshots of the 12 candidates are shown in 
Figure 4. In this figure, the ID number in- 
creases from top to bottom, and the bandpasses are 
B, V, R, i, NB816, zb, z, NB921, zr from left to right. All 
the objects are found to be detected clearly in both zb, 
z, and NB921, while some are not visible in i or NB816, 
which are expected to be at higher rcdshifts. The sky dis- 
tribution of the 12 candidates is shown in Figure 5. The 
dotted lines outline the region in which the LBG selec- 
tion is made. The area in the lower-left corner, which 
corresponds in size to the field of view of a single CCD, 
is not used in this study, since the quantum efficiency of 
the CCD in charge of this area is about two thirds those 
of the other nine CCDs and thus the data of this area are 
considerably noisy. 

We estimate a possible contamination rate for simple 
i-dropout selection. If we omit Zb — zr < 0.7 and i — zr > 
4(zb — zr) + 0.1 criteria and thus select objects based on a 
simple i-dropout criterion (i — zr > zb — zr + 1.3, similar 
to the criteria adopted in the previous studies) and non- 
detection in B,V,R, then we have ten more candidates. 
These newly selected ten objects might be LBGs at z > 6.2 
or very dusty LBGs at z < 6, but they are more likely to 
be foreground objects. If all of them are really foreground 
objects, the contamination fraction will be ~ 50%, imply- 
ing that bright single-color selected z-dropout samples for 
z ~ 6 LBGs are highly contaminated. An accurate evalu- 
ation of contamination needs spectroscopic observation of 
these objects. It should be, however, emphasized that the 
two criteria we omit in this simulation, which set an up- 
per limit to the zb — zr color of LBG candidates, not only 
reduce a possible contamination from foreground objects 
but also provide a sharp upper limit to the redshift range 
of our LBG sample. 

Dickinson et al. (2004) found five z ~ 6 candidates in 
HST ACS data of the two GOODS fields based on i - z 
color. Among them, three are brighter than z = 25.4, 
giving a surface density of 0.009 arcmin" 2 . Stanway et 
al. (2003, 2004) made a similar i-drop survey in the 
same fields. There are in total nine LBG candidates 
with z < 25.4 in their LBG samples; the surface density 
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Fig. 4. Snapshots of 12 LBG candidates. From left to right, 
B,V,R, i,NB816,ZBi 2 iNB921,2ri images are placed. The image 
size is 6" a side. North is up and east is to the left. 

is 0.026 arcmin -2 . There is one object with z < 25.4 in 
Bouwens et al.'s (2003) data of the ACS Guaranteed Time 
Observation fields (46 arcmin 2 ), giving a surface density of 
0.02 arcmin -2 . These values differ by up to a factor of two 
from ours, 12 candidates / 767 arcmin 2 = 0.016 arcmin -2 . 
In general, however, surface density measurements depend 
on selection criteria, redshift range surveyed, and detec- 
tion completeness, which are different among the authors. 
Thus, we do not make a detailed comparison of surface 
density. 

5. Results and Discussion 

5.1. Number Density of Bright LBGs 

Since the 12 LBG candidates are distributed in a nar- 
row range of zr magnitude, zr = 24.8 - 25.4, we do not 
fit the Schcchtcr function to the data to derive the best- 
fit values of the three Schechter parameters (a, M*, and 
cf)*). Instead, we measure the normalization of the lumi- 
nosity function (num mag -1 Mpc -3 ) at the average zr 
magnitude of our sample (zr = 25.1). This average mag- 
nitude corresponds to an rest-frame far-ultraviolet (FUV: 
around 1400A) absolute magnitude of Mpuv = —21.6 mag 
at z = 5.9, the central redshift of our survey. The abso- 
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Fig. 5. Sky distribution of 12 LBG candidates. LBG selec- 
tion is made in the region outlined by the dotted line. 

lute magnitude is calculated as Mpuv ~rn — 51og(di(^)) + 
2.51og(l + z), where c?l(z) is the luminosity distance. The 
rest-frame wavelength corresponding to the central wave- 
length of the zr band changes from 149lA (at z = 5.6) to 
1367A (2 = 6.2). We assume that LBGs have flat spectra 
(/„ = const) over this wavelength range. 

The value of the luminosity function at Mpuv = —21.6 
is calculated as: 



1 



21.6) [mag -1 Mpc -,J ] = — 



12 



1 



0.6 f-f V eS (rm) 



(1) 



where m-i is the zr magnitude of the i-th candidate, 
Vcff (mi) is the effective survey volume (comoving) for the 
i-th candidate, and 0.6 is the size of the magnitude bin. 
The effective volume for a given apparent magnitude m is 
computed by V e g(m) = f z p(m, z)^jjdz, where p(m,z) is 
the probability that an LBG of magnitude m and redshift 
z is detected in our zr image and passes the selection crite- 
ria defined above when photometric errors are taken into 
account (Steidel et al. 1999), and ^ is the differential 
comoving volume for a solid angle of 767 arcmin 2 . 

We estimate p(m, z) assuming that the color distribu- 
tion of z ~ 6 LBGs is similar to that of bright z ~ 4 
LBGs detected in the SDF by Ouchi et al. (2004a) and 
that color differences among LBGs are due to different 
E(B - V) values. Ouchi et al. (2004a) found that the 
E(B — V) of z ~ 4 LBGs is distributed rather uniformly 
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over < £?(i? — V) < 0.3 (They adopted a single, young 
age for their LBGs to estimate E(B — V)). We generate a 
number of young LBGs uniformly over < E(B — V) < 0.3, 
assign them apparent zr magnitudes, and distribute them 
randomly on the original zr image. The value of p(m,z) is 
calculated as the ratio of the number of simulated LBGs 
(with given zr mag and redshift z) which are detected 
and pass the selection criteria defined above, to the total 
number of simulated LBGs. In this simulation, we assume 
that the shape of LBGs is PSF-like (FWHM = 1".0). The 
systematic difference between input magnitudes and mag- 
nitudes measured by SExtractor is found to be negligible 
in our simulation. The median FWHM value of simulated 
objects measured by SExtractor is also found to agree 
with that of the 12 LBG candidates. The probability for 
zr = 25.1, p(25.1,z), is a function of redshift with a peak 
of 0.70 at z = 5.9 and a rather steep slope on both low- 
redshift and high-redshift sides; its FWHM is about 0.6. 
Hence, the relation between apparent magnitude and ab- 
solute magnitude is sufficiently tight in our sample. 
We obtain 0(-21.6) = (2.6 ±0.7) x 10" 5 mag" 1 Mpc -3 , 



where the error indicates the Poisson error (= 1/V12)- 
Although we have not searched for the best-fit values for 
the Schechter parameters, the fact that no LBG candidate 
is found at zr < 24.8 in our data appears to favor a faint 
Mp UV . About four LBGs with zr < 24.8 are expected to 
be found if we fix a = —1.6 and Mp UV = —21.04, which are 
the values for z ~ 3 LBGs (0* is adjusted so that </>(— 21.6) 
takes the correct value). Similarly, adopting Afp UV = —20 
and —22 gives 0.7 and 9 candidates, respectively. In the 
previous studies, Bouwens et al. (2004a) have obtained 
relatively faint Mp UV = —20.26 from their i-drop sample, 
while Bunker et al. (2004) found that the M* of z - 6 
LBGs is consistent with the value observed at z ~ 3. 

Bouwens et al. (2004a) derived the luminosity function 
of i-drop galaxies at z ~ 6 from HST ACS data of the 
GOODS fields and the Hubble Ultra Deep parallel fields. 
The normalization of the bright part (Mpuv < —19.7) of 
their luminosity function has been determined primarily 
from the wide-field (> 300 arcmin 2 ) GOODS data. The 
set of the best-fit Schechter parameters they obtained is 
a = -1.15, M* uv = -20.26, and 0* = 0.00173 Mpc" 3 , 
which gives <£(-21.6) =4.3 x 10" 5 Mpc" 3 . This 0(-21.6) 
value is about 70% as high as the value obtained by us. 
The major portion of this difference can be accounted for 
by the Poisson errors in both measurements, but cosmic 
variance probably contributes to it to some degree. In 
any case, our number-density measurement must be more 
reliable, since it is based on wide-field two-color selection. 
Our survey area is larger than the sum of the two GOODS 
fields by more than factor two, thus being less sensitive 
to cosmic variance. Cosmic variance is probably higher 
for brighter LBGs, since observations of LBGs at z < 5 
have revealed that brighter objects tend to be clustered 
more strongly (e.g., Giavalisco & Dickinson 2001; Ouchi 
et al. 2004b). We obtain a rough estimate of the cosmic 
variance following Somcrvillc et al. (2004). The right 
panel of their figure 3 shows that the root variance of 
dark matter at z = 6 for our survey volume (1.1 x 10 6 
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Fig. 6. FUV luminosity functions of LBGs at z ~ 3 (filled 
triangles: Steidel et al. 1999), z ~ 4 (open triangles: Ouchi 
ct al. 2004a), z ~ 5 (open circles: Ouchi et al. 2004a), and 
z ~ 6 (stars: Bouwens et al. 2004a; filled circle: this study). 

Mpc 3 ) is 0.033 in the A = 0.7 Cold Dark Matter model. 
Assuming a one-to-one correspondence between LBGs and 
dark haloes, we find from the left panel of the same figure 
that the bias parameter of dark haloes at z = 6 with the 
same number density as of our LBGs (1.5 x 10~ 5 Mpc -3 ) 
is 7.3. The cosmic variance is thus predicted to be about 
25%(= 0.033 x 7.3). 

We plot in Figure 6 the FUV luminosity function mea- 
surements for LBGs at z ~ 3 (Steidel et al. 1999), z ~ 4 
and 5 (Ouchi et al. 2004a from shallower SDF data), and 
z ~ 6 (Bouwens et al. 2004a and this study) . The number 
density of faint z ~ 6 LBGs detected in the HUDF (Yan 
& Windhorst 2004; Bunker et al. 2004) appears to be 
consistent with Bouwens et al.'s (2004a) luminosity func- 
tion. We focus on bright (M < — 21) LBGs for which our 
measurement is available. An interesting trend is seen 
that bright LBGs are less numerous at higher redshift, al- 
though the difference between z ~ 5 and 6 may not be so 
significant. The number density of bright LBGs drops by 
a factor of about five from z ~ 3 to ~ 6. 

The observed drop of the number density could be due 
to either dust absorption or a real decrease in LBGs. 
Adelberger & Steidel (2000) found no significant correla- 
tion between apparent magnitude and E(B—V) for LBGs 
at z ~ 3. A similar result was obtained by Ouchi et al. 
(2004a) for z ~ 4 LBGs. These results imply that appar- 
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ently brighter LBGs will on average have brighter intrin- 
sic luminosities, although there is a large scatter between 
intrinsic and apparent luminosities. Ouchi ct al. (2004a) 
also found that the E{B — V) distribution does not change 
over z ~ 3 and z ~ 4. If the dust properties found in z ~ 3 
- 4 LBGs apply to z > 5 LBGs as well, then we can in- 
terpret Figure 6 as suggesting that the number density of 
intrinsically bright LBGs decreases with redshift. Another 
possibility is that the observed decline in bright LBGs is 
due to an increase in the dust absorption of bright LBGs 
with redshift; in this case, the number density of intrinsi- 
cally bright LBGs may not necessarily drop with redshift. 
However, there is so far no observation showing evolution 
in dust properties of LBGs with redshift. In §5.3, we dis- 
cuss the evolution of the FUV luminosity density of bright 
galaxies over < z < 6, assuming that the dust properties 
are unchanged over this redshift range. 

5.2. Star Formation Rates 

We calculate the star formation rates of the 12 LBG 
candidates from their FUV continuum luminosity using 
the formula given in Madau et al. (1998), which assumes 
a Salpeter IMF with 0.1M Q < M < 125M Q . In the cal- 
culation, all the candidates are assumed to be located at 
z = 5.9. The results are given in Table 2. The star forma- 
tion rates of the 12 candidates are found to be ~ 20 - 30 
Mq yr _1 , with an average of 23 M Q yr _1 . The amount 
of dust extinction for the candidates cannot be estimated 
from our data, but if it is similar to the typical value for 
z ~ 3 - 4 LBGs, E(B - V) ~ 0.15, then the extinction- 
corrected star formation rate of each candidate will be 
close to 1 x 10 2 M© yr _1 . If these candidates maintain the 
observed star formation rates for about 1 Gyr since z = 6, 
they will have stellar masses of > 3 x 10 10 M Q at z = 3. 
This value is comparable to the observed stellar masses of 
bright 2^3 LBGs (Papovich et al. 2001; Shapley et al. 
2001). 

5.3. Evolution of the FUV Luminosity Density 

We measure the completeness-corrected FUV luminos- 
ity density, pl, contributed from LBGs brighter than 
Alpuv = -21.3 as: 

PL (< -21.3) [ergs s^Hz^Mpc" 3 ] = f] ffiL (2) 

r-f U e ff(mi) 

where L(rrii) is the FUV luminosity (ergs s _1 Hz -1 ) corre- 
sponding to apparent magnitude TOj. The absolute mag- 
nitude Mpuv = —21.3 is close to Mp UV of z ~ 3 — 4 LBGs. 
We find p L (< -21.3) = (2.8 ± 0.8) x 10 24 ergs s" 1 Hz" 1 
Mpc~ 3 , which is equivalent to a star formation rate den- 
sity of (3.5 ± 1.0) x 10~ 4 M o yr" 1 Mpc" 3 . If we integrate 
the luminosity function of Bouwens et al. (2004a) over 
Mpuv < -21.3, we obtain 6.3 x 10" 4 A/ Q yr" 1 Mpc" 3 . 
Stanway et al. (2003, 2004) estimated the star formation 
rate of z ~ 6 LBGs brighter than z = 25.6 in each of the 
two GOODS fields. They found (6.7±2.7) x 10~ 4 M Q yr" 1 
Mpc~ 3 for the GOODS-South and (5.4 ±2.2) x 10~ 4 Af Q 
yr -1 Mpc" 3 for the GOODS-North. These values are 



higher than ours, mainly reflecting the higher number den- 
sity of bright LBGs selected in the GOODS fields than in 
the SDF as discussed in §5.1. 

As will be discussed below, our sample shows that the 
FUV luminosity density of bright (M FUV < -21.3) LBGs 
at z ~ 6 is 11 times lower than that at z ~ 3. This ratio is 
reduced to 6 - 8 if Stanway et al.'s (2003, 2004) estimates 
are adopted. In the same GOODS fields, Giavalisco et al. 
(2004) found a factor of 3.5 decrease from z ~ 3 to z ~ 6 
for ^850 < 26 LBGs. The smaller decreasing factor found 
by Giavalisco et al. (2004) is probably due to the fainter 
limiting magnitude of their sample, since the evolution 
of FUV luminosity density tends to be milder for fainter 
galaxies (see below). 

We examine how the luminosity densities for bright 
galaxies and all galaxies evolve over < z < 6, and whether 
or not there is a remarkable difference in their evolution- 
ary behavior. Figure 7 shows the observed evolution of 
FUV luminosity density over < z < 6. The large symbols 
indicate the luminosity density contributed from bright 
(Mpuv < —21.3) galaxies, pl{< —21.3), while the small 
symbols are for the luminosity density from galaxies down 
to a faint limit of M F uv = —19.0, pl(< —19.0). We use 
Pl(< —19.0) as a substitute for the total FUV luminosity 
density in the universe, assuming that pl(< —19.0) traces 
the relative change in the total luminosity density accept- 
ably, although the contribution from galaxies fainter than 
— 19 mag is uncertain and may be considerably large. The 
filled circle indicates our measurement, while the other 
kinds of symbols represent previous measurements taken 
from the literature. We use Bouwens et al.'s (2004a) lu- 
minosity function to calculate pl(< —19.0) at z = 6. The 
vertical axis on the right-hand side of Figure 7 shows cor- 
responding star formation rate density. Dust extinction 
correction has not been applied to any measurements. The 
star formation rate density calculated from pl(< —19.0) 
multiplied with factor 10 is roughly consistent with re- 
cent measurements of the extinction-corrected total star 
formation rate density (Hopkins 2004). 

We find that pl(< —21.3) evolves drastically with red- 
shift. It increases by an order of magnitude from z ~ 6 
to rsj 3, has a rather sharp peak at z ~ 3, and then drops 
by (more than) three orders of magnitude to the present 
epoch. The evolution of pl(< —19.0), on the other hand, 
is much milder. Over 1 < z < 6, pl(< —19.0) is roughly 
constant within a factor of about five, and then it de- 
creases by ten times from z ~ 1 to the present epoch. As 
a result, Ph(< — 21.3)//jl(< —19.0) reaches a maximum 
around z = 3. 

The remarkably different evolution of Pl(< —21.3) and 
Pl{< —19.0) probably reflects a change in Mp UV with 
time. It has been found that Afp UV fades by 3 magni- 
tude from z = 3 (~ —21) to z = (~ —18) (Arnouts et 
al. 2005). A faint Mp UV is suggested for z ~ 6 LBGs 
since no candidate brighter than = 24.8 is detected. 
Since Mpuv = —21.3 is comparable to or brighter than 
Mp UV at any redshift, the luminosity function has always 
a sharp, exponential shape around Mpuv = —21.3. Thus, 
Pl (< —21.3) is much more sensitive to the change in Mp UV 
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Fig. 7. Evolution of the FUV luminosity density. The 
large and small symbols indicate the luminosity density con- 
tributed from galaxies brighter than JWprjV = —21.3 and 
Mfuv = —19.0, respectively. The filled circle with error bars 
show our measurement, and the other symbols indicate data 
taken from the literature: cross: Bouwens et al. (2004a), 
open circles: Ouchi et al. (2004a), open triangles: Steidel et 
al. (1999), open squares: Arnouts et al. (2005; HDF), filled 
triangles : Arnouts et al. (2005; GALEX), stars: Sullivan 
et al. (2000). The two filled triangles at the bottom imply 
that the actual measurements are lower than 1 X 10 22 ergs 
s- 1 Hz" 1 Mpc" 3 . 

than p h (< -19.0). 

An absolute magnitude of Mpuv = —21.3 corresponds 
to a relatively high star formation rate of 18M© yr _1 (or 
a dust-corrected value of 80M© yr _1 if E(B — V) = 0.15 is 
assumed). Figure 7 shows that such actively star forming 
galaxies become abundant in a short cosmic time from z = 
6 to 3 (about 1.2 Gyr) and then disappeared by the present 
time. Their contribution to the cosmic star formation was 
the highest at z ~ 3. 

The drastic evolution of actively star forming galaxies 
(i.e., bright LBGs) found here may be physically inter- 
preted as a combination of two competing processes in a 
Cold Dark Matter universe: the build-up of massive dark 
haloes with time and the decrease in the star formation 
rate per unit dark matter mass (i.e., specific star forma- 
tion rate) in dark haloes with time. Hernquist & Springcl 
(2003) and Springcl & Hernquist (2003) have made essen- 
tially the same discussion to explain the observed cosmic 
star formation history. 



Clustering analyses of z ~ 3 — 5 LBGs have shown that 
dark haloes hosting bright z ~ 3 - 5 LBGs have a total 
mass of Mtot ~ 1O 12 M (e.g., Giavalisco & Dickinson 2001; 
Moustakas & Somerville 2002; Ouchi et al. 2004b). Let 
us assume that bright z ~ 6 LBGs also reside in ~ 10 12 M Q 
dark haloes. According to the standard Cold Dark Matter 
model with fi = 0.3, Ao = 0.7, h = 0.7, and ag = 0.9, 
the fraction of matter in the universe confined in dark 
haloes with M tot > 10 12 M Q increases by an order of mag- 
nitude from z — 6 to 3. This is similar to the observed 
evolution rate in the luminosity density of bright LBGs. 
Note that the model predicts the fraction of matter in 
Mtot > 10 12 M Q haloes to increase by an order of magni- 
tude from z = 3 to 0; the increase per unit cosmic time is 
thus more gradual at z < 3 than at 3 < z < 6. 

The decrease with time in the specific star formation 
rate may also be qualitatively in accord with Cold Dark 
Matter models. Dark haloes formed later have a lower 
mass density, since the average mass density within dark 
haloes formed at a given cosmic time is proportional to 
the average mass density of the universe at that time. It 
is reasonable to expect that specific star formation rate 
correlates positively with the mean mass density within 
dark haloes. Hernquist & Springcl (2003) and Springcl 
& Hernquist (2003) have shown that gas cooling becomes 
insufficient in dark haloes formed at lower redshifts be- 
cause of the decline of the mass density, and thus the spe- 
cific star formation rate in dark haloes declines with time. 
They found that in their numerical simulations the cosmic 
star formation rate has a peak at z ~ 5 — 6. There may 
be, of course, other mechanisms which suppress star for- 
mation in haloes formed recently, such as the truncation 
of star formation in galaxies infalling into galaxy groups 
and clusters which become abundant with time, and the 
decrease in galaxy merging involving star formation. 

If we adopt this interpretation, the observed evolution 
of actively star forming galaxies implies that the spe- 
cific star formation rate of massive (M to t ~ 10 12 M Q ) dark 
haloes is roughly unchanged over 3 ^ z ^ 6, while it drops 
drastically with time at z < 3 for some reason, resulting 
in the star formation activity of the whole massive haloes 
being the highest around z = 3. The mild evolution of the 
total luminosity density over z ~ 3 and 6 looks reason- 
able, since the dark-matter mass summed over all haloes 
increases by only a factor of two or so from z = 6 to 3. On 
the other hand, the mild evolution at z ^ 3 may be due 
to a gradual decrease, for some reason, in the specific star 
formation rate in less massive haloes. 

Finally, we point out an interesting resemblance of the 
evolution of bright LBGs (i.e., FUV-selected actively star- 
forming galaxies brighter than —21.3 mag) to that of lu- 
minous dusty star-forming galaxies (submm galaxies) and 
QSOs. Chapman et al. (2003, 2004, 2005) found that 
the space density of luminous dusty star-forming galax- 
ies increased by about an order of magnitude from z ~ 4 
to z ~ 2 and then decreased by ~ 10 3 from z ~ 2 to the 
present time, having a sharp peak at z ~ 2.2. Chapman 
et al. (2003) have pointed out that luminous dusty star- 
forming galaxies coexist with the peak activity of QSOs. 
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Recent wide-field QSO surveys have measured accurately 
the luminosity function of bright QSOs at z ^ 6, finding 
that the space density of bright (Mpuv < —26.7) QSOs 
increases by more than an order of magnitude from z ~ 6 
to ~ 3, shows a maximum around z = 2.5, and then drops 
by about 10 3 by the present epoch (Croom et al. 2004; 
Fan et al. 2001, 2004). 

These three populations appear to be related with ac- 
tive star formation or massive galaxy formation. Dusty 
star-forming galaxies, detected in submillimeter or ra- 
dio observations, have star formation rates of as high as 
- 10 3 M© yr" 1 (e.g., Chapman et al. 2003, 2004). It is 
widely believed that QSOs are the result of gas accretion 
onto a supermassive black hole at the center of a galaxy, 
and hence are likely be linked to active star formation in 
the galaxy and, possibly, to bulge formation (Franceschini 
et al. 1999). Nearby luminous QSOs are hosted by mas- 
sive bulge-dominated galaxies (e.g., Floyd et al. 2004). 
It has also been found that the host galaxies of luminous 
QSOs at 1 < z < 2 are consistent with massive galaxies 
following passive evolution (e.g., Falomo et al. 2004). 

At z ~ 2 - 3, the number density of bright LBGs, ~ 10~ 4 
Mpc~ 3 , is higher than those of dusty star- forming galaxies 
(~ 10" 5 Mpc" 3 ) and QSOs (~ 10" 7 Mpc" 3 ). This may 
imply either that the lifetime of bright LBGs, i.e., the du- 
ration of star formation of > 20 — 30A/ Q yr -1 , is longer 
than the lifetimes of dusty star-forming phase and QSO 
phase, or that bright LBGs reside in less massive (thus 
more numerous) dark haloes than the other two popula- 
tions, or both. A weak clustering of z ~ 3 LBGs compared 
with the clustering of QSOs (Croom et al. 2005) and 
dusty star-forming galaxies (Blain et al. 2004) at similar 
redshifts could imply that LBGs reside in the least mas- 
sive (but still as massive as ~ 10 12 M Q ) dark haloes. The 
relatively milder evolution in number density at z > 2 for 
LBGs than for dusty star-forming galaxies may also sug- 
gest smaller dark halo masses of LBGs. It is also interest- 
ing to note that the peak redshift for bright LBGs (z ~ 3) 
is slightly higher than those for dusty star-forming galax- 
ies and QSOs (z ~ 2 - 2.5); this might reflect hierarchical 
evolution of haloes that less massive haloes are formed 
earlier. 

The physical relationship between the three populations 
is not yet clear, but the resemblance of their evolution sug- 
gests that z ~ 3 is a remarkable era in the cosmic history 
when massive galaxies were being intensively formed. The 
cosmic time of z ~ 3 also falls in the suggested formation 
epoch of massive early- type galaxies seen at < z < 2 (e.g., 
Bower, Lucey, & Ellis 1992; Kodama et al. 1998; Stanford 
et al. 1998; van Dokkum & Stanford 2003; Cimatti et al. 
2004; McCarthy et al. 2004). 

6. Conclusions 

We showed that LBGs at z ~ 6 are isolated well from red 
galaxies at z ~ 1 - 2 and Galactic cool stars in the i — zrvs 
z~b — z~r plane. We applied this two-color selection to wide- 
field (767 arcmin 2 ) multicolor data of the Subaru Deep 
Field, selecting 12 candidates for z ~ 5.9 ±0.3 LBGs down 



to zr = 25.4. This is the most reliable sample of bright 
z ~ 6 LBGs, which is robust against both contamination 
by foreground objects and cosmic variance. 

The average star formation rate of the 12 candidates is 
23 Mq yr -1 if dust extinction is not corrected. We cal- 
culated the normalization of the rest-frame far-ultraviolet 
(FUV: ~ 1400A) luminosity function of z ~ 6 LBGs at 
M FUV = -21.6 to be 0(-21.6) = (2.6±0.7) x 10~ 5 mag" 1 
Mpc -3 . The FUV luminosity density contributed from 
LBGs brighter than M FUV = -21.3 is (2.8 ±0.8) x 10 24 
ergs s _1 Hz -1 Mpc -3 after correction for detection com- 
pleteness, which corresponds to a star formation rate den- 
sity of (3.5 ± 1.0) x 10- 4 Af o yr- 1 Mpc~ 3 . 

We examined the evolution of the FUV luminosity den- 
sity, or equivalently the star formation rate density, over 
< z < 6 by combining our measurement with those taken 
from the literature. We found that the FUV luminos- 
ity density for bright FUV galaxies with Afpuv < —21.3, 
i.e., actively star- forming galaxies, increases by an order 
of magnitude from z ~ 6 to ~ 3 and takes a rather sharp 
peak at z ~ 3, followed by 10 3 times decrease toward z = 0. 
On the other hand, the evolution of the 'total' luminosity 
density was found to be much milder. As a result, the 
contribution to the cosmic star formation from actively 
star-forming galaxies is the highest at z ~ 3. The dras- 
tic evolution of actively star-forming galaxies can be ex- 
plained by a combination of the build-up of massive dark 
haloes with time and the decrease in the specific star for- 
mation rate in dark haloes with time. We also pointed 
out a resemblance between the evolution of bright LBGs, 
luminous dusty star-forming galaxies, and bright QSOs. 
The redshift of z ~ 3 seems to be a remarkable era in 
the cosmic history in terms of the formation of massive 
galaxies. 
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Table 2. Photometric properties of LBG candidates. 



ID 


a 


(J2000.0) 


6 (J2000.0) 


ZR 


zb - ZR 


i - 


" ZR 


NB816 


z 


NB921 


SFR (a) 


1 


13 


23 


46.55 


+27 


27 


3.9 


24.84 


0.12 


2 


.41 


27.04 


25.96 


25.88 


29.2 


2 


13 


25 


27.81 


+27 


28 


58.9 


24.92 


0.65 


2 


.79 


27.04 


25.42 


25.37 


27.1 


3 


13 


24 


39.15 


+27 


32 


13.0 


24.94 


0.47 


2 


.50 


27.04 


26.06 


26.42 


26.6 


4 


13 


23 


49.72 


+27 


45 


27.4 


25.02 


-0.14 


1 


.47 


26.10 


26.57 


26.07 


24.7 


5 


13 


24 


16.91 


+27 


13 


48.3 


25.04 


0.06 


2 


.28 


27.04 


26.01 


26.16 


24.2 


6 


13 


25 


23.15 


+27 


42 


51.1 


25.14 


0.26 


2 


.66 


27.04 


26.24 


25.76 


22.1 


7 


13 


24 


9.43 


+27 


18 


30.4 


25.15 


0.36 


1 


.85 


25.75 


25.55 


25.42 


21.9 


8 


13 


23 


39.77 


+27 


31 


37.4 


25.16 


-0.01 


1 


.89 


27.04 


25.49 


25.68 


21.7 


9 


13 


24 


27.13 


+27 


16 


40.1 


25.23 


0.48 


2 


.22 


27.04 


26.13 


25.95 


20.4 


10 


13 


23 


55.87 


+27 


32 


9.3 


25.26 


0.14 


1 


.61 


26.05 


25.73 


25.94 


19.8 


11 


13 


24 


40.01 


+27 


43 


32.3 


25.36 


0.32 


1 


.93 


26.95 


26.19 


25.87 


18.1 


12 


13 


25 


7.62 


+27 


45 


7.4 


25.36 


0.43 


2 


.49 


27.04 


25.99 


25.83 


18.1 



